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Abstract: We consider the PDWHEG generator for a H/W/Z boson plus one jet, augmented 
with the recently proposed MiNLO method for the choice of scales and the inclusion of 
Sudakov form factors. Within this framework, the generator covers all the transverse- 
momentum region of the H/W/Z boson, i.e. no generation cuts are needed to obtain a finite 
result. By construction, the generator achieves NLO accuracy for distributions involving a 
finite (and relatively large) transverse momentum of the boson. We examine the conditions 
under which also the totally inclusive distributions (e.g. the boson rapidity distribution) 
achieve NLO accuracy. We find that a minimal modification of the MiNLO prescription is 
sufficient to achieve such accuracy. We thus construct a NLO generator for H/W/Z boson 
plus one jet production such that it smoothly merges into a NLO single boson production 
in the small transverse- momentum region. We notice that, by simply reweighting the boson 
rapidity distribution to NNLO predictions, we achieve a NNLO accurate generator matched 
to a shower. The approach applies to all production processes involving a colorless massive 
system plus one jet. We discuss how it may be extended to general processes. 
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1. Introduction 

In recent times, next-to-leading order parton shower (NLO+PS) matching techniques have 
been developed and realised as practical simulation tools [1, 2, 3]. The programs and 
methods underlying them have matured to the point where they are routinely used in LHC 
data analysis. 

NLO+PS methods can also be applied to processes involving associated jet production 
at leading order. In the case of Higgs boson production, for example, there exist POWHEG 
simulations for inclusive production of a Higgs boson (the H generator from now on) [4], 
Higgs boson production with an associated jet (HJ), and also for the case of two associated 
jets (HJJ) [5]. These simulations overlap in their population of phase space. However, 
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the relative accuracies of each one in the various regions is complementary. Specifically, 
the H generator yields NLO accurate inclusive Higgs boson distributions, describing the 
radiation of a single jet with leading-order (LO) accuracy and the radiation of more than 
one jet with the accuracy of the shower to which it is interfaced (the collinear, leading-log, 
approximation). On the other hand, the HJ generator gives NLO precision for observables 
inclusive in the production of a Higgs boson plus one jet; it has only LO quality if we require 
two jets, with further radiation generated in the parton shower approximation. The HJ 
generator cannot be used to describe inclusive jet cross sections, since the NLO calculation 
on which it is based requires the presence of at least one jet. Similar considerations hold 
for the HJJ generator. 

It is natural to ask whether one can merge the H, HJ and HJJ simulations, in such a 
way that all classes of observables have NLO accuracy in the end results. That is to say, 
one would ideally like to have a unified NLO+PS simulation, or simulation output, yielding 
inclusive Higgs boson distributions accurate at the NLO level, and, at the same time, 
distributions involving the Higgs boson plus one jet accurate at order af. 

Several collaborations have addressed this merging problem [6, 7, 8, 9, 10, 11, 12, 13]. 
They all typically separate the output of each component simulation (the equivalent of 
H, HJ or HJJ) according to the jet multiplicity of the events it produces, discarding those 
having a multiplicity for which the generator does not possess the relevant NLO corrections. 
Having processed the output of each simulation in this way, the event samples are joined 
to give an inclusive sample. In a nutshell, each generator can therefore be regarded as 
contributing a single exclusive jet bin to the final inclusive sample, the magnitude of each 
bin being predominantly determined by the jet resolution scale used in performing the 
merging, the so-called merging scale. The merging scale is an unphysical parameter, and the 
dependence on it is, rightly, well studied in the implementations of the merging algorithms 
cited above, with different approaches invoking different means to mitigate the dependence 
on it. Relatedly, we point out that in all of the practical applications of NLO merging 
to date, the jet algorithms employed do not correspond to those currently used in LHC 
experimental analyses, albeit for good theoretical reasons: in particular, to lessen the 
dependence on the merging scale through the resummation of spurious large logarithmic 
corrections associated with it. 

In all cases the choice of the merging scale poses a dilemma: if the merging scale is too 
low, the sample is dominated by the higher-multiplicity generators, while if it is chosen too 
large, one loses their benefit, since one is forced to describe relatively hard jets only with 
tree-level accuracy, or, worse still, with the parton-shower approximation. 

To further clarify the situation consider again Higgs boson production. As stated 
above, the H generator yields 0(a|) accuracy for inclusive quantities, i.e. quantities that 
are integrated over the transverse momentum of the Higgs boson. Consider, however, the 
same integrated cross section up to a certain transverse momentum cut p^ ut . Assuming 
that the Sudakov form factor provided by the NLO+PS implementation is next-to-leading 
log (NLL) accurate, missing next-to-next-to-leading log (NNLL) corrections to this quantity 
(corresponding to a factor ~ expfagL] ~ 1 + a^L + where L = log(M H /p!j, ut ), and M H is 
the Higgs boson mass) reduce the accuracy of the distribution, depending upon the size of 
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a s L. If we assume a s L 2 ~ 1, as in the vicinity of the Sudakov region, these neglected NNLL 
terms are of relative order a^' 5 with respect to the Born term. Thus, the NLO quality (that 
requires neglected terms to be of relative order a 2 or higher) is lost. Thus, although the 
inclusive distributions are NLO accurate, the cut cross section may well not be if the cut is 
too low. Of course, if the cut is high enough the NLO precision is restored, since NLO+PS 
simulations describe the high transverse-momentum tail with tree-level accuracy and the 
fully integrated spectrum at the NLO level. 

Turning now to look at the HJ generator from the same perspective, if we consider the 
cross section with a transverse momentum cut p T > plj, ut , we know that large logarithmic 
corrections of the form a s L 2 relative to the Born term arise in the cross section, hence, 
if a s L 2 ~ 1, the predictions completely fail. Some of the new NLO merging techniques 
attempt to resum these logarithmic corrections to all orders, in general though this resum- 
mation is, with the notable exception of ref. [11], limited to NLL contributions and with 
the neglect of the NNLL terms, as in the H case, leading to errors at the level of ~ «g' 5 
terms relative to the leading-order cross section. 

In order to avoid this problem, the authors of ref. [10], conservatively, use merging 
scales which are not too low. The authors of ref. [13] proposes to tackle the problem by 
forcing unitarity on the approach, using suitable subtractions in order to restore NLO 
accuracy. In ref. [11] the accuracy of the resummation is improved in order to reach the 
required precision for matching. 

In the present work we study the possibility of building a NLO generator that does 
not lose NLO accuracy when the radiated parton is integrated over. In other words, we 
want to achieve the goals of merging without actually doing any merging at all. 

In a recent paper [14], some of us have proposed a prescription (MiNLO: Multi-scale 
Improved NLO) for the choice of renormalization and factorization scales and the inclusion 
of Sudakov form factors in NLO calculations, such that the higher-multiplicity calculations 
seem to merge well with the lower-multiplicity ones. The numerical results of ref. [14] are 
quite striking, in that they seem to suggest that with this method, for example, the HJJ 
simulation yields a good description also for observables that are best computed with the H J 
or H programs. In ref. [14], however, no arguments are given to justify this behaviour, and 
no attempt is made to quantify the formal accuracy that a generator which is NLO accurate 
in the description of (n + l)-jet observables has for those sensitive to m jet topologies, with 
m < n. 

In the present work, we address precisely these issues. We consider BJ generators, 
where B denotes either the Higgs (H), the W (W) or the Z boson (Z) and ask the following 
questions: 

• How precise is the BJ generator, when improved with the MiNLO prescription, in 
describing inclusive boson production observables? 

• Is it possible to modify the MiNLO procedure in such a way that we achieve NLO 
accuracy also for these observables? 

We will show that: 



- 3 - 



• The inclusive boson observables are described by the BJ+MiNLO programs at relative 
order a s with respect to the Born cross section. However, they do not reach NLO 
accuracy, since they also include ambiguous contributions of relative order a| 5 , rather 
than «g. 

• The second question has a positive answer. In the main body of the paper we will 
clarify what needs to be done in order to gain genuine NLO precision for inclusive 
observables. The modifications needed for reaching this goal are simple and we have 
implemented them in our current MiNLO-POWHEG BOX HJ, WJ and ZJ generators. 

• As a final point, we notice that by reweighting the new BJ-MiNLO generators in such a 
way that the inclusive boson rapidity and decay kinematics are matched to the next- 
to-next-to- leading order (NNLO) result, one obtains a NNLO calculation matched to 
a parton shower simulation, i.e. a NNLO+PS generator. 

The event generators built in this way can be further enhanced, so as to yield an improved 
description also for the third hardest emission, as proposed in ref. [7]. To this end, it 
is enough to simply feed the events obtained with the BJ-MiNLO programs through the 
relevant BJJ simulation, with the latter attaching the third hardest parton to the events 
with LL+LO accuracy. 

The paper is organized as follows. In Section 2 we discuss the precision of the original 
BJ+MiNLO generator and the modifications that are necessary to attain NLO accuracy in 
the description of inclusive boson distributions. In Section 3 we give an alternative, more 
detailed proof. In Section 4 we present results obtained with our revised MiNLO approach 
and discuss issues related to the scale variation as an estimator of the theoretical error on 
the predictions. In Section 5 we show how we can improve the generators that we have 
constructed in order to achieve NNLO accuracy. In Section 6 we discuss briefly how the 
present findings may be generalized to more complex processes. Finally we present our 
conclusions in Section 7. 

2. Accuracy of the BJ+MiNLO generators 
2.1 Preliminaries 

We want to determine the accuracy of the BJ+MiNLO generators when integrated over the 
boson transverse momentum at fixed rapidity. We will illustrate explicitly the case of Higgs 
boson plus jet production, but our findings hold trivially for WJ and ZJ as well, and for the 
production of a heavy colour- neutral system accompanied by one jet. 

In order to avoid confusion, we will refer to the accuracy of inclusive results, integrated 
over the boson transverse momentum, as the (0) accuracy, and will refer instead to the 
accuracy in the H + j inclusive cross section as the (1) accuracy. Thus, for example, by 
LOW and NLO(°) we mean leading (i.e. «g) and next-to-leading (i.e. ccf) order accuracy 
in the Higgs boson inclusive rapidity distribution, and with LO^ and NLO^ we mean 
leading (i.e. af) and next-to-leading (i.e. a|) order accuracy in the Higgs boson plus one 
jet distributions. 
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We focus upon the HJ generator with the MiNLD prescription. As illustrated in [14], this 
is obtained by modifying the POWHEG B function with the inclusion of the Sudakov form 
factor and with the use of appropriate scales for the couplings, according to the formula 



B = a 2 s {Ml) a s (q 2 T ) A 2 g (M H , q T ) B (l - 2A« (M H , q T j) + V + J d$ rad R 



(2.1) 



where we have stripped away three powers of a s from the Born (B), the virtual (V) and 
the real (R) contribution, factorizing them in front of eq. (2.1). We will comment later on 
the scale at which the remaining power of a s in R, V and is evaluated. 



A„ is the gluon Sudakov form factor 



A g (Q, q T ) = exp 



rQ 



f[A(a s ^))lo g f + B(a s ( q 2 )) 



and 



A g (Q,q T ) = 1 + A« (Q,q T ) + 0(a 2 s ) 



(2.2) 



(2.3) 



is the expansion of A g in powers of a s - In the previous equations, q T stands for the Higgs 
boson transverse momentum in the underlying Born kinematics, and Q 2 its virtuality. 
The functions A and B have a perturbative expansion in terms of constant coefficients 



(2.4) 
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In MiNLO, only the coefficients A\, A2 and B\ are used. They are known for the quark [15] 
and for the gluon [16] and are given by 
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The O (q s ) expansion of the Sudakov form factor in eq. (2.3) is given by 

1 



A« (Q,q T )=a s 



,2 1 

T 



2 ^log 2 | + ^log Q2 



(2.5) 
(2.6) 

(2.7) 
(2.8) 

(2.9) 



In ref. [14], we give a particular prescription for the choice of the renormalization scale in 
the power of a s accompanying V, R and A g *\ and for the explicit renormalization and 
factorization scales present in V and R, that we do not need to specify at this moment. 
Now we will assume, at variance with [14] (for reasons that will become clear later), that 
the renormalization scale for the power of a s accompanying V, R and A^ is q T . 
Our argument is based upon the following considerations: 
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1. Equation (2.1) yields the exact NLO^ 1 -* cross section if expanded in terms of the 
strong coupling (two powers of them evaluated at the scale M H , and the remaining 
ones at the scale q T ), and of the parton densities evaluated at the scale q T . Thus, 
when integrated in the whole phase space, keeping fixed only the Higgs boson q T and 
rapidity, and expanded up to order a|, eq. (2.1), must include, for small q T , all the 
singular parts of the NLO^ result for the H + j cross section. 

2. The singular parts of the H + j cross section at NLO^ can also be obtained using 
the NNLL resummation formula for the Higgs boson transverse momentum. 

3. The NNLL formula, plus the regular part of the LO^ H + j cross section, when 
integrated over all values of q T , is NLO^ ) accurate. This is due to the fact that, up 
to order a| (M^) , the NNLL formula includes terms of the following form: 



\ Qt J + \ Qt / + 

and thus, when adding to it the regular terms, we achieve NLO^ accuracy. 

4. By comparing eq. (2.1) and the NNLL formula, we can thus see what is missing in 
eq. (2.1) in order to achieve NLO^ accuracy. 

Before we continue, however, we should remember that, in the MiNLO formula, q T stands 
for the Higgs boson transverse momentum in the underlying Born kinematics, rather than 
the true Higgs boson transverse momentum. We ignore this difference, assuming that a 
similar formula also holds for the underlying Born transverse momentum. Later on we will 
modify the MiNLO prescription to achieve full consistency. 

We adopt the following NNLL formula for the Higgs boson transverse-momentum 
distribution at fixed rapidity 

xex P 5(Q,^ T )} + i? / , (2.11) 

where f a /A and f b / B are the parton distribution functions for partons a and b in hadrons 
A and B, respectively, xa and xb denoting the partons' momentum fractions at the point 
of annihilation (y = ^log^-). The convolution operator, 0, and Sudakov exponent, S, 
are defined 

(f®g)(x)= f 1 -f(z)g(-), exp<S (Q, q T ) = Ag (Q, q T ) , (2.12) 

j x z \ z / 

with the coefficient functions Cij having the following perturbative expansion 

oo 

Cij (a s ,*) = <M(1 - z) + Y j a n s C^(z) . (2.13) 

n=l 
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Lastly, we have used Rf to label the non-singular part of the cross section. The renor- 
malization and factorization scales in eq. (2.11) are set to q T , as indicated explicitly in the 
argument of the convolution. 

Transverse-momentum resummation is usually performed in impact-parameter space. 
A formulation of resummation in q T space for vector-boson production of the form (2.11) is 
given in ref. [17]. It is however integrated in rapidity. Furthermore it is not of NNLL, and 
not even of NLL order in the usual sense [18]. In Appendix A we show that it is however 
possible to derive a complete NNLL resummation formula in q T space starting from the 
6-space formula. It will become clear in the following that the difference of this formula 
from the one of ref. [17] is in fact irrelevant for our proof. 

In the following, together with the explicit values of A\, A2 and Bi, we need the 
expression for the B2 coefficient too. For Higgs boson production, it was computed in 
ref. [19]. Including the modification needed to go from the impact-parameter space to the 
transverse- momentum expression (see ref. [17] and also Appendix A) it assumes the form 1 



B: 



fl(H) 



1 



2 " 2tt 2 

+ 4C 3 K) 2 , (2.14) 
while for Drell-Yan production it was computed in [20, 21] and is given by 

< Y) -h [(£ - h- *) * + - £ + 1*) ™ + (s - a) H 

+ 4CsK) 2 - (2-15) 
2.2 The NLO(°) accuracy and BJ+MiNLD 

We now examine in more details under which conditions eq. (2.11) achieves NLO(°) accu- 
racy. By integrating it in dg 2 , we get 



da 
dy 



^0 [C ga (8) fa/A (xa, Q) x [C gb ® f b/B ] (x B , Q) + J dq\ Rf , (2.16) 



where we have neglected the lower bound of the integration of the first term, since, at small 
<7tj the integrand is strongly suppressed by the Sudakov exponent, exp5 (Q,q T ). We thus 
see that, in order to reach NLO*- ^ accuracy, the C g % and C g j functions should be accurate 
at order a s (i.e. they should include the term) and Rf should be LO*- 1 -* accurate. 

Notice that the form of eq. (2.11) (i.e. the fact that it is given as a total derivative) 
is such that the NLO(°) accuracy is maintained by construction, independently of the 
particular form of the Sudakov form factor, as long as one includes the terms. 

We now want to show that even if we take the derivative in eq. (2.11), and discard 
terms of higher order in a s , the NLO(°) accuracy is maintained. We consider eq. (2.11) 
after the derivative is taken, for q T . In order to maintain its relative NLO^ 1 ^ accuracy, we 



1 Notice that, because of the expression of the Sudakov form factor in eq. (2.12), our A and B coefficients 
are divided by 2 with respect to those of refs. [19, 20, 21]. 
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need to include the C\- terms, the A\ and A2 terms, and the B\ and B2 terms of eq. (2.4). 
After the derivative is taken, we get terms of the following form 

a \ [a s , a?, af, aj, a s L, a 2 s L, a^L, ajL] expS {Q,q T ) , (2.17) 
Qt 

where L = log Q 2 /q 2 and the a s L and a 2 L terms arise from the A term in the derivative 
of the Sudakov exponent. Some terms of order a s and a 2 also arise in this way, from 
the B term. Others arise from the derivative of the parton distribution functions (pdfs). 
Terms of order a| arise, for example, from two terms together with the B2 term in the 
derivative of the Sudakov exponent. All powers of a s in the square bracket of eq. (2.17), 
and the relevant pdfs, are evaluated at q T . 

If we do not drop any higher-order terms from eq. (2.17), its integral is still NLO^°) 
accurate. This is, of course, the case, since the formula can be written back as an exact 
derivative, its integral is given by eq. (2.16), and the terms are included. We now want 
to show that even if we drop the higher-order terms, the NLO^ ** accuracy is not spoiled. 

We estimate the size of each contribution using the formula 

^ da 2 O 2 r, ra+l 

/ ^ log™\<(<£) exp5(Q, gT )« [a s (Q 2 )] n ~ . (2.18) 
J y\2 q T q T 

This formula is a consequence of the fact that the dominant Sudakov singularities carry 
two logarithms for each power of a s , and thus each logarithm counts as We give 

however a detailed derivation of this formula, including the effect of the running coupling, 
in Appendix C. 

Among all terms of order higher than a 2 in the square bracket of eq. (2.17), the 

3— - 

dominant one is a^L. This term gives a contribution of relative order [a s (Q 2 )]' 2 = 
a 2 (Q 2 )- Thus, all terms of order a| and higher can be dropped without spoiling the 
NLO(°) accuracy. Dropping these terms, we get essentially the full singular part of the 
MiNLO HJ formula, except that the original MiNLD formula does not have the B2 term in 
S. In fact, if we expand the Sudakov factor up to 0(a s ) and combine it with the content 
of the squared parenthesis, we get the full singular part of the HJ cross section. Since also 
the MiNLO HJ formula has the same property, the two must agree. We also remark that the 
choice of the scale in the power of a s entering V, and R (that we have taken equal 
to q T ) is essential for our argument to work. For example, if we choose instead a scale 
equal to Q, the largest difference arises in the terms of order a 2 L in eq. (2.17), where they 
would give an O (aigL 2 ) variation. This yields a contribution of order O («g' 5 (Q 2 )) upon 
integration. On the other hand, we will show in the following this kind of contributions 
are of the same size as the effect of the B2 term. Thus, if the latter is not included, this 
scale choice remains ambiguous. 

We now investigate what is the loss of precision due to the lack of the B2 term in 
the MiNLO formula. In order to do this, we drop the B2 term in the Sudakov exponent 
in eq. (2.11). The resulting formula still satisfies eq. (2.16), and thus is NLO^ ) accurate. 
When taking the derivative, however, we will get all the terms in the square bracket of 
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eq. (2.17), except for the term 

a \a 2 s {q 2 )B 2 expS(Q,q T ). (2.19) 
Qt 

This term is instead present in the MiNLO result, that agrees with eq. (2.17) up to the 
terms of order a 2 in the square bracket, and differs from it only by subleading terms in the 
Sudakov exponent. We conclude that, if we take away the contribution of eq. (2.19) from 
the MiNLO result, we recover the NLO(°) accuracy. Thus, the MiNLO formula, as is, violates 
the NLO(°) accuracy by this term, which yields, according to eq. (2.18), a contribution of 
relative order O (a^ 5 (Q 2 )) upon integration. 

2.3 Summary 

Summarizing our findings: 

• the original BJ+MiNLO generator of [14], is less than NLO^ accurate, in that it 
includes incorrect terms of relative order O (og' 5 (Q 2 )). 

• In order to achieve NLO^ accuracy for the BJ+MiNLO generator, we must 

- include the B 2 term in the Sudakov form factor of eq. (2.2); 

— take the scale of the power of a s entering V, and R equal to q T . 

In this section we have not considered the possibility of varying the factorization and 
renormalization scales in the MiNLO improved B functions. It is indeed possible to vary 
them in such a way that the NLO^ accuracy for the BJ+MiNLO generator is maintained. 
In Appendix B we show in detail how this can be done. 

3. Secondary proof of NLO accuracy 

In this section we describe an independent derivation, corroborating that of Sect. 2, that 
MiNLO heavy boson plus jet production computations, with minor modifications, yield pre- 
dictions NLO accurate in the description of both inclusive and jet-associated observables. 
The alterations to the original MiNLO algorithm required here are the same as those dis- 
cussed previously in Sect. 2. Principally this amounts to using the transverse momentum 
of the heavy boson, as opposed to the scale obtained on clustering the Born kinematics 
with the & T -jet algorithm, as an input to the pdfs, strong coupling and Sudakov form 
factor, together with the inclusion of the NNLL B 2 coefficient in the latter. The analysis 
essentially comprises of four steps: 

1. Using existing expressions for their singular parts, we write the NLO cross section for 
heavy boson plus jet production, differential in the boson's rapidity, y, and transverse 
momentum, q T , as the sum of a q T — > divergent part, da s , and a finite remainder, 
do>, with scales set according to the slightly revised MiNLO conventions. 
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2. The MiNLO Sudakov form factor and coupling constant weights are applied to the 
resulting expression and the related subtraction terms, necessary to maintain NLO 
accuracy in the presence of these weights, are inserted. 

3. By neglecting terms which, on integration over q T , give rise to unenhanced terms 
O («g) with respect to the lowest order heavy boson production process, we are able 
to cast the singular part of the MiNLO cross section in the form of a total derivative. 

4. The q T integration of the singular cross section may be performed at NLO accuracy 
and combined with the leading part of the remainder cross section, do>, whereupon 
the sum can be identified as the NLO differential cross section. 

In the following we give details regarding these steps, pointing out any deviations with 
respect to the original MiNLO algorithm as they enter. The mathematical description of the 
algorithm, as given here, reflects, precisely, that of the implementation whose results we 
present in Sect. 4. Lastly, we further stress that although the formulae we use explicitly 
refer only to the heavy bosons' rapidity distributions, the derivations hold, with a simple 
modification, also for the case of distributions involving the heavy-boson decay products. 
This generalization is discussed in Sect. 3.4. 

3.1 NLO q T spectra with MiNLO scale choices 

To begin proving the NLO accuracy of the MiNLO predictions for inclusive observables, we 
require an expression for the cross section differential in the transverse momentum and 
rapidity of the heavy boson, employing the renormalization and factorization scales set out 
in ref. [14]. Contrary to the original MiNLO prescription, however, we now insist that the 
transverse momentum of the massive boson be used as the factorization scale, moreover, 
we now also wish that this scale, q T , be used as the argument in evaluating the extra a s 
factor accompanying the NLO corrections (ag LO ). 

The asymptotic q T — > limit of W and Z transverse momentum spectra was given in 
ref. [22] to next-to-leading order accuracy. The NLO Higgs boson transverse momentum 
spectrum has also been computed in ref. [23]. In the latter work the authors derived the 
limit directly from their NLO computation and found agreement with the expression ob- 
tained by taking the analogous limit of the 6-space resummation formula of ref. [24]. For 
W, Z and Higgs boson production processes the results are quoted for arbitrary renormal- 
ization (/%) and factorization scales (/i F ). 

Using these results one can directly write down NLO expressions for the small trans- 
verse momentum limit of the heavy-boson q T spectra, by simply replacing all instances 
of /% and fi F (explicit and implicit) by K R Q and K F q T , respectively, in the asymptotic 
expressions of refs. [22, 23]. K F and K R simply denote constant rescaling factors used to 
perform scale uncertainty estimates, thus they are nominally set equal to one, nevertheless, 
we retain them explicitly in order to better clarify the nature of such variations. Having 
made the aforesaid replacements, one can trivially substitute the additional power of a s in 
the NLO terms by ag LO = a s (K^Q 2 ) simply, with no further changes of any kind. 
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We thus write the q T spectrum of a heavy boson resulting from the annihilation of 
partons i and j , from beam particles A and B respectively, as the sum of a q T — > singular 
contribution, da s , and a finite regular remainder term, do>, 



do- 



da s 



da 7 



dq-^dy dq-^dy dq-^dy ' 
where the singular part takes, precisely, the form 

2 2r- 1 



(3-1) 



da s 
dq^dy 



N 



(l)"EE( 



r=l s=0 



(3.2) 



with L = log (Q 2 /q^) ■ The r D s coefficients are directly obtained from the rn C n coefficients 
given in Appendix A. 2 of ref. [22] and Appendix C of ref. [23], following the replacements 
described above. Since the resulting expressions are lengthy but trivial to derive we refrain 
from quoting them explicitly here. Lastly, in the normalization of eq. (3.2) the factor of a™ 
accounts of there being n powers of the strong coupling associated with the leading order 
production process (n = for W and Z production, n = 2 for Higgs boson production) 
while A comprises constant factors such that the product of A(a s /27r) n with the pdfs 
gives the LO cross section for H/W/Z production in the i, j channel, differential in y. 

Equation (3.1) describes, precisely, the cross section returned by the MiNLO programs 
prior to the introduction of Sudakov form factor and coupling constant reweightings. In 
particular, since the cross section we write here is intended to refer to that coded in the 
programs, there are no unknown O (a"" 1 " 3 ) terms or beyond omitted anywhere, throughout 
it contains explicitly only contributions proportional to and a™ +2 , with all strong 

coupling and pdfs factors utilizing the scale settings already elaborated on. 



3.2 Differential cross section in the MiNLO programs 

Following the MiNLO procedure we must now multiply the NLO cross section by a Sudakov 
form factor and, to maintain NLO accuracy, simultaneously subtract a term corresponding 
to the O («s +2 ) contribution implicit in the product of the leading order terms with the 
form factor. In so doing, as was with the choice of factorisation scale, here we wish to 
change the lower scale entering the Sudakov form factor, from that obtained on the first 
clustering of the underlying Born kinematics with the /c T -jet algorithm, to simply q T . 

In our programs we have used the following Sudakov form factor, describing the evolu- 
tion of a quark / gluon line from scale Q to q T , generalised to account for a renormalization 
scale rescaling factor K R 



log A (Q, K R q T ) 



Q 2 



wherein 



S , 



-.2 



A(a s (q*))lo g 9l + B(a s (q>)) 



and 



S , 



(3.3) 



(3.4) 
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with the Ai and Bj coefficients related to the A and B; of Sect. 2 as follows 



i 2 = A + A & log i^R 



B 2 



Bi + A x log if 2 



+ ^n& 2 bg K R 2 + A. log Kl+ l -A l b log 2 X| . 



(3.5) 
(3.6) 
(3.7) 

(3.8) 



Here, in the definition of B 2 , we have used B { ^ to denote either B^ H ^ or B\^ as 
defined in eqs. (2.14) and (2.15) respectively. The derivation of the scale dependence of 
these coefficients is given in Appendix B. 

For the case of NLO Higgs or vector boson plus jet production processes, the MiNLO 
algorithm will multiply the full cross section by the product of two such Sudakov form 
factors. Introducing the abbreviation, L = L — logK 2 , the first order expansion of this 
aggregate form factor is given by 



? 9(DY) 



A 2 (Q,K R q T ) = l-a s 



Ai L 2 + 2 Bi L 



(3.9) 



To compensate the spurious O (a™ + ) contribution which this will generate, when it mul- 
tiplies the leading order part of the heavy boson plus jet cross section (distinguished by 
the subscript LO), a corresponding counterterm is added. 

In addition to the Sudakov form factor, the MiNLO procedure, for Higgs, W and Z 
plus jet processes, prescribes that the whole cross section be multiplied by an overall factor 
comprised of a ratio of coupling constants 



a s (K 2 q 2 ) /a s {kIQ 2 ) = l + a s b L + O {a 2 s L 2 ) . 



(3.10) 



As for the factorisation scale and the low scale for the Sudakov form factor, here we have 
revised the original MiNLO algorithm by taking the scale in a s in the numerator to be K R q T . 
To apply this overall weight factor and maintain the integrity of the NLO cross section, a 
counterterm is introduced in the MiNLO cross section, balancing the erroneous contribution 
which it elicits at O {a™ +2 ). 

Taking into account the Sudakov form factor, the coupling constant ratios, and their 
respective counterterms, the MiNLO cross section, da M , reads as follows 



dq^dy 



A 2 (Q,K R q T ) 



da 



+ 



a s (K^g 2 ) 
a s (KIQ 2 ) 
da 



dg^dy dq^dy 



a. 



LO 



IiL 2 + 2B ± L - b (L + log 



(3.11) 



All modifications explained in this cross section have been defined precisely, with no am- 
biguities of any kind, higher order or otherwise, that is to say, eq. (3.11) reflects, exactly, 
the implementation of the differential cross section in the MiNLO programs for Higgs, W 
and Z plus jet production. 
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Combining eqs. (3.1), (3.2) and (3.11), using only basic algebraic manipulations, we 
may write the cross section from the MiNLO programs in such a way as to isolate the leading 
q T — > behaviour. Specifically, we further develop our da M formula as 



da M 
dq^dy 



dcr. 



da, 



dq^dy 



dq^dy ' 



(3.12) 



with the leading q T — > part given by 
d<7 Ms 



dq^dy 



^ T r=l s=0 



S r-'-'s 



(3.13) 



and subleading part given by eq. (3.11) save for the replacements da M —> da MT on the 
left- and dcr — > do> on the right-hand side. In eq. (3.13) n powers of a s are evaluated at a 
scale K R Q and the remainder at K R q T . Assuming the same shorthand as ref. [22] for the 
pdfs, fi = fi/A (xa,^f), the r E s coefficients are found to be 



i Ei — 2 Ax fi fj , 

iEq = 2B\ fi fj + [P ik 

2E X = A 2 fi fj + 2A 1 

2 E = B 2 fi fj + 25i 



(3.14) 



fk] fj + [Pjk ® fk] fi , 
/* fj + j} , 



°ik 



ik 



fk fj + 



□ (2) 
ik 



fk fj + boln^ [P lk 



+ 



cT®fk [Pn®fi] + 



C 



ik 



kl 



fi fj ~ h 



a 



ik 



fk] fj 

fk fj j] , 



with Cjfc" defined, exactly, as 



C ik" i z ) = C \k (*) + ^ b log Kl S tk 5 (1 - z) - log P tk (z) 



•<(i) 



Ai log 2 Kl + B 1 log Kl 6 ik 6 (1 - z) . 



(3.15) 



The C\l ] terms are NLO corrections to the coefficient functions in the conventional 6-space 
resummation formula, while Pi k (z) denotes the leading order splitting function for parton k 

(2) 

branching to a parton i with momentum fraction z, P- k representing its relative O (a s ) 
corrections. The and P ik functions here are equal to those in refs. [22, 23] divided by 
2tt, similarly, our pjj^ functions are equal to those in refs. [22, 23] divided by (27r) 2 . 2 

In determining eq. (3.12) from eq. (3.11) we have at no point employed Taylor expan- 
sions, renormalization group or DGLAP equations; eq. (3.12) follows exactly from eq. (3.11) 
without neglect of terms at O (a™ +3 ) or beyond. Thus eq. (3.12) corresponds precisely, 
without ambiguities of any kind, higher order or otherwise, to the differential cross section 
returned by the MiNLO programs. 



2 In simplifying the r E s coefficients, in particular for collecting the C^}" terms, the following trivial 
identities are useful: fi = SikS (1 — z) ® fk and P ik (g) fk = SikS (1 — z) ® Pm ® fi ■ 
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3.3 Integrating the MiNLO q T spectrum with NLO accuracy 

In this subsection, through judicious omission of higher order contributions, together with 
application of the DGLAP and coupling constant evolution equations, we write the leading 
q T — > part of the MiNLO cross section, da MS , in the form of a derivative with respect 
to q T . To this end we use the following result derived from the integral in Appendix C: 
terms in the differential cross section O (a^ +r L s / 'q^) , with r > \ {s + 5), multiplied by the 
Sudakov form factor, yield contributions O (ag +2 ) or higher on integration over the low q T 
domain, i.e. such terms do not affect inclusive observables at the NLO level, nor, for r > 3, 
do they affect NLO accuracy of jet-associated H/W/Z production observables. 

Neglecting contributions not greater than O (a™ +3 L/g 2 ) times the Sudakov form fac- 
tor, one can replace the sum of terms proportional to Ai and Bi in da MS by the q T derivative 
of the Sudakov form factor, multiplied by a product of the coefficient functions, C" k , defined 

C'{ k z) = 5 tk 5 (l-z) + a^°c£" (z) . (3.16) 

Further, dropping terms O (as +3 ) inside the summation of eq. (3.13), one may substitute 

a s {Klql) (p ik + aT° (p£> + b P lk log |f ) ) ® f k -+ q\ M , (3 .17) 

with the implicit scale in the pdfs being K F q T as usual. 3 The same replacement without 
the ag LO terms on the left-hand side holds, at the same level of approximation, for the two 
remaining terms in 2-E0 containing the leading order splitting functions. Lastly, discarding 
additional terms of the same order as those omitted in eq. (3.17), one may also replace the 
remaining 60 factor in the 2-E0 coefficient as 

60^ -4 ^ (3^) 
ai dq$ 

with the scale for a s in the denominator and the derivative being K R q T . 

With the help of these substitutions the leading q T — > part of the MiNLO cross section, 
neglecting terms no greater than ~ A 2 as +3 L/q^, can finally be written as 



A M - n (S)" 4 { a2 «2- K Al K ^'] } ■ < 3 - 19 » 



where the n powers of a s explicit in the prefactor retain K R Q as their argument, with 
K n q T holding for the rest. Integrating eq. (3.20) over q T yields 



da 



MS 



(hj = N ( g)™ A 2 (Q, K R Q) [C'{ k ® f k ] [C% ft] , (3.20) 

with the coupling constants in the C'[ k having K R Q as their argument and the scale in the 
pdfs set now, here, to K F Q. The Sudakov form factor here contains no large logarithms 



3 Among the O (e*s +3 ) terms forgone in this replacement are contributions proportional to 
a" +i l°g 2 K^/Kp, which we neglect in the understanding that the ratio K R /K F will be of order 1. In 
fact, in our phenomenology section 4, the ratio is confined to the interval [|,2] . 
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since the integral range in the exponent is limited to high values, [K^Q 2 , Q 2 ~\, thus we 
may reliably expand eq. (3.20) in powers of a s neglecting NNLO sized contributions to 
give simply 



da? 



dy 



N 



2^ 



n ( 

{fifj 



+ «S 



c$? ® fk 



fj + «s 



c 



(1) 



jk 



h 



fi} 



(3.21) 



wherein all pdf and a s factors are evaluated at the common scale Q. The disappearance of 
all factors of Aj, Bi in this expression follows from the cancellation of the NLO expansion of 
the Sudakov factor, A 2 (Q,K R Q), with an equal and opposite contribution in the product 
of the C'( k coefficient functions. Similarly, the vanishing of all K R and K F factors attributes 
to the explicit logK 2 /F terms being absorbed in the pdfs and strong coupling constant 
factors, renormalizing the scales in these factors, in the leading order contribution, from 
K r /fQ to Q. The scales of the pdfs and coupling constants in the remaining (NLO) terms 
are freely changeable at this level of accuracy. 

To proceed we must explain the nature of the NLO component of the coefficient func- 
tions, . These are given, essentially as a matter of definition, by half of the sum of the 
virtual and real, q T — > singular, components of the NLO cross section for heavy boson 
plus jet production, with the real-radiation phase space below q T = Q having been inte- 
grated over. This calculation is carried out explicitly in ref. [25], for the case of Higgs boson 
production, and in Sect. 2 of ref. [26] for the qq channel in vector boson production. The 
coefficients obtained agree precisely with those in refs. [23, 25], i.e. with those employed 
in the expressions for the NLO (7 T -spectra from which we started our analysis in Sect. 3.1. 
Additionally, at the end of Appendix C in ref. [27] the integral of the W and Z q T spectra 
up to some arbitrary q T is shown, wherein one can readily identify the NLO Cy functions 
as the combined virtual and integrated real, q T — > singular, contributions to the weak 
boson cross section. 4 

Finally we assert that the integral of the remaining part of the q T spectrum, da MT , is 
equal to that of the original q T — > regular component, do>, (with the coupling constant 
and pdf scales of order Q), up to formally relative order a 2 terms. This statement follows 
from the fact that the regular parts are suppressed by a factor q 2 /Q 2 relative to their 
diverging counterparts in da s . Thus, in the region q T < Qa s {Q 2 ) surrounding the Sudakov 
peak, the regular piece is effectively a relative a 2 (Q 2 ) contribution and hence entirely 
negligible from the point of view of our NLO analysis. Just outside this region the same 
principle applies, namely, that provided the proportionate contribution made by the regular 
piece is small with respect to the finite ones, the impact that the Sudakov form factor 
makes on them is of higher order significance. This argument holds a fortiori given that 
the Sudakov form factor is introduced with compensating terms, multiplying all of the 
Born cross section, cancelling its effects at O (a s ), effectively deferring them to yet lower 
regions in the q T spectrum, i.e. to regions where the q 2 /Q 2 suppression of the regular cross 
section have already rendered its contribution small. Finally, by way of a more quantitative 



4 Up to an irrelevant notational difference of a factor of two, the coefficient functions in ref. [27] are 
identical to those of ref. [25] . 
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ratification of this claim we note the following integral, 

= a s + (ajjj) , 

(3.22) 

performed assuming the coupling to be fixed at a scale O (Q). This integral has the form of 
those which would be encountered in evaluating the da MT , showing it to be equivalent to 
the analogous do> integral. 5 Although we have resorted to a fixed coupling approximation 
here, experience suggests this gives results compatible with other approaches, moreover, for 
NLO accuracy we only require the do> and da MT integrations be the same up to O (ai?) 
corrections rather than the observed O (af ) . We therefore conclude that the contributions 
from do> and da MT may be safely taken as being equal. 

With the correspondence between the Cy coefficient functions and the NLO fixed 
order computations clear, as well as that between the q T integral of do> and da MT , we 
have thus proved that their combination, viz. the integrated / inclusive MiNLO cross section, 
is equivalent to that of a conventional fixed order computation at NLO accuracy. 

3.4 Extension to vector-boson decay products 

While our discussion has only referred explicitly to the rapidity distributions of the W, Z 
and Higgs bosons, the proof that these distributions are NLO accurate straightforwardly 
extends to the case where the cross section is further differential in the kinematics of 
their decay products. For the case of a scalar Higgs boson this is trivial, with the decay 
essentially decoupled from the production stage, occurring isotropically in its rest frame. 

In the case of vector-boson production we are not aware of explicit results in the lit- 
erature for the q T — > limit of the NLO cross sections, including the dependence on the 
momenta of the decay products. However, we may infer such expressions using ref. [28], 
which clarifies how to incorporate variables parametrizing the momenta of such particles 
in the 6-space resummation formalism. This is achieved by modifying the resummation 
formula of ref. [24] via the insertion of the unit normalised LO partonic cross section, 
differential in the decay variables, at the front of the expression (outside the 6-space in- 
tegration), together with the inclusion of a hard function, H^, multiplying the coefficient 
functions Cij. The hard function partly encodes the NLO corrections to the leading order 
distribution of the decay products' momenta. Clearly, with these localised changes, if one 
carries out a perturbative expansion of this 6-space formula one will re-obtain the asymp- 
totic q T — > cross sections of refs. [22, 23], modified by an overall factor of the leading 
order cross section differential in the decay variables and with the O (a s ) terms in the 
coefficient functions modified to include NLO corrections to that distribution. Thus in- 
cluding dependencies on the details of the vector-boson decay leaves our proof unmodified 
up to inclusion of an overall angular factor describing the decay kinematics and, in general, 
an implicit dependence on these variables in the O (a s ) part of the coefficient functions. 

Taking a different tack, for the case that the decay variables are the Collins-Soper 
angles [29], this extension of the resummation formalism was already well studied in 

5 The same integral without the exponential or log 2 Q 2 /q 2 terms. 
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[Q 2 dq 2 q 2 / 



Q 2 



2 a s I 1 + a s A log" — 2~ | exp 



) exp - 2 L 



— A a s log — 2 



transverse-momentum space [17] and impact-parameter space [27, 30]. In the former ar- 
ticle, in transverse momentum space, one can readily see the angular dependence of the 
resummed cross section present as a global angular prefactor. In fact, in refs. [30, 27, 17] 
the only modification of the 'standard' resummation formulae is just that: the coefficient 
functions are not modified with respect to those used for resummation of vector-boson 
production processes where the decay is not considered. We attribute this to the fact that 
the virtual QCD corrections, simply act to rescale the hadronic tensor for the production 
vertex by an overall factor, i.e. the dependence on the Collins-Soper decay angles (defined 
in the decaying boson's rest frame) is the same in the virtual and Born contributions to 
the cross section. Moreover, only those parts of the real cross section which are singular 
when q T — > are integrated over in combining with the virtual corrections to give the 
NLO Cij terms, as one would expect, these parts therefore have an angular dependence 
for the decay products which is the same as that of the Born term. 6 Consequently, the 
only modification needed to the asymptotic q T spectrum formula in ref. [22], to account 
for these angular variables, is the inclusion of an overall factor describing the dependence 
at the relative O («g) / Born level. 

4. Implementation and plots 

In the following we discuss results obtained from an implementation of our revised MiNLD 
method in the PDWHEG BOX, for the HJ [4, 5], WJ and ZJ generators [31, 32]. For the 
purpose of this study, we will use throughout the pdf set MSTW2008NLO [33]. Any other 
popular set [34, 35] can be used equivalently, and we are not interested here to explore pdf 
dependencies. For the Higgs boson case, we use throughout M H = 125 GeV and in the H 
generator we set the hfact parameter to 100 [36]. 

We have advocated that the MiNLO prescription, updated as described in the previous 
section, can achieve NLO^ accuracy, that is to say, it can describe inclusive boson distri- 
butions at NLO accuracy. We thus begin with the most inclusive quantity, i.e. the total 
cross section. We would like to see if the total cross section obtained by integrating the 
BJ-MiNLO formulae is compatible with the NLO total cross section obtained with standard 
NLO calculations. Notice that we expect agreement only up to terms of higher order in 
a s , since the BJ-MiNLO results include terms of higher order, and also since the meaning of 
the scale choice is different in the two approaches. For similar reasons, we do not expect 
the scale variation bands to be exactly the same in the two approaches. 

In tab. 1, we list the results for the total cross sections obtained with the HJ-MiNLO and 
the H programs, both at full NLO level and at leading order, for different scale combinations. 
The scale variation in the HJ-MiNLO result is obtained by multiplying the factorization scale 
and each of the several renormalization scales that appear in the procedure by the scale 
factors K F and K R , respectively. The Sudakov form factor is also changed according to 
the prescription of Appendix B. The maximum and minimum of the cross-sections are 
highlighted. In the case of MiNLO, the leading-order result is obtained by keeping only the 

6 This is elaborated on in Appendices B and E of ref. [27]. 
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Born term in POWHEG (i.e. by setting the bornonly flag to 1), and by downgrading the 
Sudakov form factor to pure NLL accuracy, i.e. we set B2 to zero. 

In the MiNLO case, the central value is chosen according to the procedure discussed 
earlier, with more than one renormalization scale for each phase space point. In the H 
fixed order calculation, we choose as central renormalization and factorization scales the 
boson mass. From the table, it is clear that the standard NLO result and the integrated 



Higgs boson production total cross sections in pb at the LHC, 8 TeV 


K R ,K F 


1,1 


1,2 


2,1 






1 1 
2 5 2 


2,2 


HJ-MiNLO NLO 


13.33(3) 


13.49(3) 


11.70(2) 


13.03(3) 


16.53(7) 


16.45(8) 


11.86(2) 


H NLO 


13.23(1) 


13.28(1) 


11.17(1) 


13.14(1) 


15.91(2) 


15.83(2) 


11.22(1) 


HJ-MiNLO LO 


8.282(7) 


8.400(7) 


5.880(5) 


7.864(6) 


18.28(2) 


17.11(2) 


5.982(5) 


H LO 


5.741(5) 


5.758(5) 


4.734(4) 


5.644(5) 


7.117(6) 


6.996(6) 


4.748(4) 



Table 1: Total cross section for Higgs boson production at the 8 TeV LHC, obtained with the 
HJ-MiNLO and the H programs, both at full NLO level and at leading order, for different scales 
combinations. The maximum and minimum are highlighted. 



HJ-MiNLO one are fairly consistent, both at the NLO and at the LO level. At the NLO 
level, the renormalization-scale variation dominates the uncertainty band, and it turns out 
to be very similar for the HJ-MiNLO and H results, with the first one being slightly shifted 
upwards. The central values are even closer. Notice that the factorization scale variation 
is wider for the HJ-MiNLO result, a fact that we will comment on later. 

At leading order the HJ-MiNLO central result exceeds the fixed order one by almost 
50%. We again see that the renormalization scale variation dominates the uncertainties. 
The scale variation, however, is quite larger than that of the fixed order result. 

For W~ production we have considered both the LHC at 8 TeV configuration (tab. 2) 
and the Tevatron at 1.96 TeV (tab. 3). Here we notice that the WJ-MiNLO NLO result 



W — > e v production total cross sections in nb at the LHC, 8 TeV 


K R ,K F 


1,1 


1,2 


2,1 






1 1 

2 ' 2 


2,2 


WJ-MiNLO NLO 


4.35(1) 


4.65(1) 


4.031(7) 


3.818(8) 


4.84(2) 


4.62(2) 


4.462(8) 


W NLO 


4.612(8) 


4.738(8) 


4.552(8) 


4.425(7) 


4.687(8) 


4.530(8) 


4.703(8) 


WJ-MiNLO LO 


3.182(1) 


3.862(1) 


2.713(1) 


2.4531(1) 


5.006(2) 


3.792(2) 


3.305(1) 


W LO 


4.002(6) 


4.379(7) 


3.999(6) 


3.566(6) 


3.999(6) 


3.566(6) 


4.379(7) 



Table 2: Total cross section for W~ — > e~v production at the 8 TeV LHC, obtained with the 
WJ-MiNLO and the W programs, both at full NLO level and at leading order, for different scales 
combinations. The maximum and minimum are highlighted. 



has a much wider scale-variation band than the fixed-order one. In both cases, the band 
is larger by about a factor of 3. The central value is lower in both cases by about 4-5%. 
In the leading order case, the WJ-MiNLO scale band is more than twice as large as the fixed 
order one at the LHC. At the Tevatron, the scale variation for the W LO result is clearly 
too small, the NLO result being incompatible with it. On the other hand, for both LHC 
and Tevatron predictions, if only symmetric scale variations are considered (i.e. the last 
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W — > e v production total cross sections in nb at the Tevatron 


K R ,K F 


1,1 


1,2 


2,1 


M 




i i 

2' 2 


2,2 


WJ-MiNLO NLO 


1.204(2) 


1.225(2) 


1.155(1) 


1.177(2) 


1.240(4) 


1.288(4) 


1.186(1) 


W NLO 


1.253(1) 


1.264(1) 


1.232(1) 


1.246(1) 


1.281(1) 


1.273(1) 


1.241(1) 


WJ-MiNLO LO 


0.9633(3) 


0.9988(3) 


0.8334(3) 


0.8962(3) 


1.5971(6) 


1.4696(5) 


0.8677(3) 


W LO 


1.127(1) 


1.130(1) 


1.127(1) 


1.116(1) 


1.127(1) 


1.116(1) 


1.130(1) 



Table 3: Total cross section for W~ — > e~v production at the 1.96 TeV Tevatron, obtained with 
the WJ-MiNLO and the W programs, both at full NLO level and at leading order, for different scales 
combinations. The maximum and minimum are highlighted. 



two columns of the tables), the WJ-MiNLO NLO band contains the fixed order one, and is 
comparable to it. 

We notice now, that comparing full independent scale variation in the WJ-MiNLO and 
in the W approaches does not seem to be totally fair. In fact, in the W case, there is no 
renormalization scale dependence at LO, while there is such a dependence in WJ-MiNLO. 
Since the W and WJ-MiNLO cross sections are related in a known way, we can track the 
effects of scale variation in both formulae, at least in the LO case. We write, schematically, 
the leading order W production cross section as 

a (K F ) = [C (KlQ 2 ) ® cr ] , (4.1) 

where with the convolution sign we schematically represent the integral with the luminosity 
£. We have explicitly indicated the factorization scale dependence in the formula. Here Q 
stands for the reference scale, that, in this case, is the W mass. In order to make contact 
with the LO MiNLO formula, we rewrite eq. (4.1) as 

r Q2 d 

a (K F ) = dq 2 —^ e W S(Q,q) [C {K 2 q 2 ) ® a ] , (4.2) 
Jo a Q 

where 

S(Q, q) = - J ^a s (ft 2 ) (A! log ^ + B^j , (4.3) 

that, being the integral of an exact differential, yields eq. (4.1) up to non-perturbative 
effects, that we here neglect, related to the low-end of the integration. Taking explicitly 
the derivative, we get 



a(K F ) = J ^expS(Q,q)L s (Kjq 2 ) [P ® £ (K 2 F q 2 )] 

+ a s {q 2 ) (A! log ^+B^jC {K 2 q 2 ) j a 

rQ 2 ( \ f p ( r 

= / -\expS(Q,q) \a s (K 2 q 2 ) P ® £ (K 2 q 2 ) 
Jo 1 I L 



+ [a s {q 2 ) - a s {K 2 q 2 )] log ^ + B^j £ (K 2 q 2 ) ® a | . (4.4) 
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We have made use of the schematic equation 
[P <g> C (Kjq 2 )] + 



o 2 

A, log \ + B l 
q 2 



C {K 2 q 2 ) = [P®£ {K 2 q 2 ) 



(4.5) 



where P stand for the regularized (i.e. including the plus prescriptions) Altarelli-Parisi 
splitting functions, while P are the unregularized ones. In the last case, we assume that 
the integration in the convolution has the correct kinematic cutoff. 

On the other hand, the MiNLO formula for the total cross section has the form 

rQ 2 do 2 r r - i 1 

a (K R , K ¥ ) = / -\a s {Klq 2 ) P®C {K 2 q 2 ) exp5(Q, q) a + C {K 2 q 2 ) ® R f \ , 
Jo Q 1 1 J } 

(4.6) 

where the last term corresponds to the Rf finite contribution. We thus see that an inde- 
pendent scale variation in the W formula corresponds at least in part to a symmetric scale 
variation in the MiNLO formula. It is thus not surprising that the MiNLO independent scale 
variation is so much larger than the W one also at NLO. If we limit ourselves to consider only 
symmetric scale variations, the MiNLO and the W results are more consistent, although the 
W scale variation band is extremely small. The full NNLO result for the W cross section, 
on the other hand, lies outside the W uncertainty band in this case, which suggests that, 
after all, the MiNLO scale band is not unreasonable. 

In tables 4, 5 and 6 we show the total cross section for Z — > e + e~ production at the 
Tevatron and at the LHC at 8 and 14 TeV. As in the W case, we see that the full 7-point 



Z — > e+e production total cross sections in nb at the Tevatron 


K R ,K F 


1,1 


1,2 


2,1 


b| 


|,1 


1 1 

2 ' 2 


2,2 


ZJ-MiNLO NLO 


0.2407(4) 


0.2434(3) 


0.2321(3) 


0.2378(3) 


0.2459(8) 


0.2552(7) 


0.2359(3) 


Z NLO 


0.2498(3) 


0.2513(3) 


0.2455(3) 


0.2491(3) 


0.2552(3) 


0.2543(3) 


0.2467(3) 


ZJ-MiNLO LO 


0.19508(7) 


0.19883(7) 


0.16948(6) 


0.18572(6) 


0.3213(1) 


0.3029(1) 


0.17336(6) 


Z LO 


0.2052(2) 


0.2041(2) 


0.2052(2) 


0.2052(2) 


0.2052(2) 


0.2052(2) 


0.2041(2) 



Table 4: Total cross section for Z~ -4- e + e~ production, obtained with the ZJ-MiNLO and the 
Z programs, both at full NLO level and at leading order, for different scales combinations. The 
maximum and minimum are highlighted. 



Z — > e + e production total cross sections in nb at the LHC, 8 TeV 


K R ,K F 


1,1 


1,2 


2,1 


bl 


It 


1 i 

2 ' 2 


2,2 


ZJ-MiNLO NLO 


1.044(2) 


1.111(2) 


0.974(1) 


0.933(1) 


1.157(3) 


1.116(3) 


1.067(1) 


Z NLO 


1.107(1) 


1.135(1) 


1.093(1) 


1.066(1) 


1.125(1) 


1.091(1) 


1.126(1) 


ZJ-MiNLO LO 


0.7751(3) 


0.9260(3) 


0.6632(2) 


0.6111(2) 


1.2114(4) 


0.9405(3) 


0.7950(3) 


Z LO 


0.9597(9) 


1.0443(1) 


0.9604(9) 


0.8647(9) 


0.9604(9) 


0.8647(9) 


1.0443(1) 



Table 5: Total cross section for Z~ — > e + e~ production, obtained with the ZJ-MiNLO and the 
Z programs, both at full NLO level and at leading order, for different scales combinations. The 
maximum and minimum are highlighted. 



scale variation yields a much wider band in the ZJ-MiNLO result with respect to the Z one. 

After having shown that the total cross sections obtained with the BJ-MiNLO generator 
are in good agreement with the standard NLO cross sections, we would like to show that 
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Z — > e+e production total cross sections in nb at the 14 TeV LHC 


K R ,K F 


1, 1 


1,2 


2,1 






1 l 

2 ' 2 


2,2 


ZJ-MiNLO NLO 


1.916(5) 


2.065(6) 


1.776(2) 


1.662(3) 


2.18(1) 


2.022(6) 


1.987(3) 


Z NLO 


2.039(3) 


2.100(3) 


2.015(2) 


1.938(2) 


2.068(3) 


1.984(2) 


2.092(3) 


ZJ-MiNLO LO 


1.3827(5) 


1.7322(6) 


1.1806(4) 


1.0348(3) 


2.1280(7) 


1.5677(5) 


1.4831(5) 


Z LO 


1.793(2) 


2.014(2) 


1.793(2) 


1.555(2) 


1.793(2) 


1.555(2) 


2.014(2) 



Table 6: Total cross section for Z~ — > e + e~ production, obtained with the ZJ-MiNLO and the 
Z programs, both at full NLO level and at leading order, for different scales combinations. The 
maximum and minimum are highlighted. 
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Figure 1: Comparison between the H+PYTHIA result and the HJ-MiNLO+PYTHIA result for the 
Higgs-boson rapidity distribution at the LHC at 8 TeV. The left plot shows the 7-point scale- 
variation band for the H generator, while the right plot shows the HJ-MiNLO 7-point band. 

also the rapidity distributions are in good agreement. We thus show in fig. 1 the rapidity 
distribution of the Higgs boson at the 8 TeV LHC, computed with the H and with the HJ- 
MiNLO generators, both interfaced to PYTHIA 6 [37] for shower. We have used the Perugia-0 
tune of PYTHIA (that is to say, PYTUNE(320) ). Hadronization, underlying event and multi- 
parton collisions were turned off. The two plots show the scale-variation band for each 
generator. The band is obtained as the upper and lower envelope of the results obtained 
by setting the scale factor parameters (K R ,K F ) to (1, 1), (1,2), (2, 1), (1, 1), (~, 1), (1, 1) 
and (2,2). We see considerable agreement between the two approaches, with the scale- 
variation band of the HJ-MiNLO result being slightly larger. 

In figs. 2 and 3 we show the Higgs transverse momentum distributions. We begin 
by noticing that the central values of the H and HJ-MiNLO generators are in very good 
agreement. This is not a surprise, since in the H generator, the parameter hfact, that 
separates the real cross section contribution into the sum of a singular and a finite one, 
was set to the value M H /1.2, motivated by the fact that this yields better agreement with 
the NNLO result. 

We notice that, for large transverse momenta, the HJ-MiNLO generator has a smaller 
scale variation band with respect to the H one. We expect this behaviour, since the HJ- 
MiNLO generator achieves NLO accuracy for one-jet inclusive distributions, while the H 
generator is only tree-level accurate. We also notice that the scale uncertainty band of 
HJ-MiNLO widens at small transverse momentum. This behaviour is also expected, since, 
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Figure 2: Comparison between the H+PYTHIA result and the HJ-MiNLO+PYTHIA result for the Higgs 
boson transverse-momentum distribution. The bands are obtained as in fig. 1. 
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Figure 3: Same as fig. 2 for a different range. 



in that direction, we approach the strong coupling regime. Observe also that the H result 
does not show a realistic scale uncertainty in the p^, < M H region. This too is understood, 
and it follows from the fact that this region is dominated by 5-type events (see refs. [36, 38] 
for a detailed explanation). 

As a last point, we see from fig. 3, that a noticeable difference in shape is present in 
the very small transverse-momentum region. This again does not come as a surprise, since 
the POWHEG-generated Sudakov form factor in the H generator differs by NNLL terms, and 
also by non-singular contributions, from the HJ-MiNLO one. Notice also that, unlike in the 
H case [38] , the scale variation in the HJ-MiNLO generator induces a change in shape of the 
transverse momentum spectrum in the Sudakov region, leading to a better understanding 
of the associated uncertainty. 

We now turn to the case of W~ production. Motivated by the discussion given for 
the total cross section case, we consider only a 3-point scale variation, i.e. K R = K F = 
{1/2,1,2} for the WJ-MiNLD generator. In fig. 4 we show the l~ rapidity distribution at 
the Tevatron computed with the W and WJ+MiNLO generators. We essentially see no shape 
difference in this distribution, therefore, as for the inclusive cross section, we find that the 
WJ+MiNLO central value is about 5% below the W one. The WJ band is slightly larger than 
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Figure 4: Comparison between the W+PYTHIA result and the WJ-MiNLO+PYTHIA result for the 
l~ rapidity distribution at the Tevatron. The bands are obtained as in fig. 1 for the W+PYTHIA 
generator, while for the WJ-MiNLO+PYTHIA generator they are obtained by taking the upper and 
lower envelope of the curves computed with K R = K F = {1/2, 1, 2}. 



the W one for central rapidities, widening towards larger rapidities. 
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Figure 5: Comparison between the W+PYTHIA result and the WJ-MiNLO+PYTHIA result for the W~ 
transverse-momentum distribution at the Tevatron. The bands are obtained as in fig. 4. 

In figs. 5 and 6 we present predictions for the W~ boson transverse-momentum spec- 
trum. In this case we find noticeable shape differences between the W and WJ+MiNLO 
distribution, especially at low . In particular, we observe that the WJ+MiNLO Sudakov 
form factor peaks at a lower value of p^ . These differences do not come as a surprise, since 
this distribution is described only at LO by the W generator, while the WJ+MiNLO descrip- 
tion is NLO accurate. We also note that the W uncertainty band is small and uniform in 
the whole p^ range. This is a feature of P0WHEG when no separation is performed between 
singular and regular contributions to the cross section. In this case, the size of the scale 
variation amounts to a factor 1 + O(ag), that is clearly too small in the moderate to large 
p^ region, where the W generator is only tree-level accurate. The error band given by the 
WJ+MiNLO generator is of an acceptable size at large transverse momenta, while it seems 
to be excessively small in the very low transverse momentum region. 

Finally, we discuss the case of Z production at the 14 TeV LHC. In fig. 7 we show the 
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Figure 6: Same as fig. 5 for a different p™ range. 




Figure 7: Comparison between the Z+PYTHIA result and the ZJ-MiNLO+PYTHIA result for the Z 
rapidity distribution at the LHC at 14 TeV. The bands are obtained as in fig. 4. 

Z boson rapidity distribution. As in the case of W production, we note that the Z J+MiNLO 
central value is lower than the Z one and that its uncertainty band is comparable at central 
rapidities, widening in the forward-backward region. We also notice a slight change of 
shape in the extreme rapidities, with the Z result remaining compatible with the Z J+MiNLD 
uncertainty band. 

In figs. 8 and 9 we show the transverse- momentum distribution of the Z boson. We 
notice the same features already observed in the W case. However we also observe now a 
considerable difference of the Z and ZJ-MiNLO distributions at large transverse momenta. 
Following arguments in ref. [39] we suggest that this difference arises due to NLO corrections 
in the ZJ-MiNLO generator which are not present in the Z one, for example those related to 
threshold logarithms. 

In conclusion, we have seen that in the case of vector-boson production [Z and W), 
the agreement of the BJ-MiNLO and the B generators is less than perfect. We also have 
shown that there is not a close correspondence between the scale variations in the two 
types of generators. We may conclude from this observation that either the scale variation 
in the MiNLO generator is excessive, and should be limited (for example by considering only 
3-point scale variations), or that the scale uncertainty obtained in the W and Z generators 
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Figure 8: Comparison between the Z+PYTHIA result and the ZJ-MiNLO+PYTHIA result for the Z 
transverse-momentum distribution at the LHC at 14 TeV. The bands are obtained as in fig. 4. 




Figure 9: Same as fig. 8 for a different range. 



underestimates the true error. We also remark that if we use the full 7-point scale variation 
for the WJ and ZJ generators, the error bands for the transverse-momentum distributions 
do not shrink as much at small transverse momenta, thus leading to better compatibility 
with the W and Z generators distributions. On the other hand, by doing so, the error on 
the total cross section might be too conservative. The problems with the application of 
our method to these generators may be related to the fact that the bulk of the Sudakov 
region is relatively low, the peak being below 5 GeV. We have ignored some problems that 
may be relevant in this region. For example, the pdf evolution switches from 5 to 4 and to 
3 flavours in this region, while we have used constant nf = 5 in the BJ-MiNLO generator. 
Similarly, inaccuracies in the pdf evolution in this region may also be relevant. While we 
have checked that our results are not affected in a relevant way by the infrared cutoffs of 
the calculation, we believe that further work will be needed to better assess the importance 
of low-p T Sudakov region. 

On the contrary, in the Higgs boson case, we find from our study a fairly good agree- 
ment between the H and HJ-MiMLO generators, both for the total cross sections and for the 
distributions. The scale variation bands are also comparable, which confirms that our re- 
sult yields an improved and more accurate description of Higgs boson production, fulfilling 
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also, in a practical, way our goals. 



5. How to build an NNLO+PS generator with MiNLO 

We argue now that a generator built according to the MiNLO prescription, improved with 
the findings of the present work, can be easily turned into an NNLO+PS generator. For 
the sake of simplicity, we consider the Higgs boson generator. We have concluded that 
such a simulation achieves C(a|) accuracy for all distributions involving at least one jet 
and 0(al) accuracy for inclusive distributions. We denote by 

d y/Hj 

the inclusive Higgs boson rapidity distribution obtained with this event generator. We 
denote by 

(r) < 5 - 2 > 

V a V J NNLO 

the inclusive Higgs boson rapidity distribution calculated with a fixed NNLO calculation. 
We claim that by reweighting the HJ-MiNLO output with the weight factor 



NNLO 



d ^/HJ 



(5.3) 



we achieve full NNLO accuracy for our generator. The proof is simple. The above ratio 
has a formal expansion 

— I 

<Wnnlo _ c 2 a s + c 3 a s + c 4 a s ^ l + c*-*^ + ^ ^ 



(da\ c 2 a| + c 3 a| + d^aj c 2 

\ d y Jej 

Notice that the numerator and denominator agree at C(a|) since the HJ-MiNLO generator 
achieves C(a|) accuracy for inclusive observables. Thus, the reweighting factor does not 
spoil the CKg accuracy of the HJ-MiNLO generator in the one-jet region. In this region, in fact, 
the dominant contributions to the HJ-MiNLO generator are of order cuf , and the reweighting 
generates extra contributions of order q|, that are beyond the nominal accuracy. On the 
other hand, the inclusive distributions are reweighted to achieve CKg accuracy, so that the 
generator has indeed a| accuracy in the whole phase space. 

Variants of these schemes are also possible. Rather than reweighting using the full 
rapidity distribution, one can split the cross section as 

d 2 <7 d 2 <7 , . , d 2 cr ri , . x , 

■ M?t) + -r-r [1 - M?t)] , (5-5) 



dg T dy dq T dy dg T dy 

where h is a smooth positive function such that h(q T ) — > 1 as q T — > and h(q T ) — > for 
q T » M H , such as 

cM a 

h(q) = „ H , (5.6) 
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with a > 1 and c a constant of order 1. One can then reweight the cross section as 

,l " 7 h(q T )W(y) + -^-[l-h(q T )] , (5.7) 



dq T dy dq T dy 

with 



dV 

IT-(y) = - vdgTdj/y NNLO (5 8) 

d cr \ 



J dq T h (q. 
dq T h (q T ) 



dq T dy) 



HJ 



In this way, for q T » M H the effect of the reweighting vanishes. 

We notice that this NNLO+PS generator would be NNLO accurate in the same sense 
in which the current MCSNLO or POWHEG type generators are NLO accurate, i.e. integrated 
quantities achieve NLO/NNLO accuracy, while LL/NLL/NNLL accuracy is achieved in the 
Sudakov region, depending upon the accuracy of the implementation of the Sudakov form 
factors. 

We postpone a phenomenological study including reweighting to a future publication. 

There are a number of more complex processes to which this procedure can be gen- 
eralized, typically processes where the heavy particle decays, or processes involving pair 
production of massive colourless objects. In this case, the NNLO reweighting must be per- 
formed as a function of more variables. One would typically use the rapidity of the heavy 
system, plus the kinematical variables describing its internal structure. For example, one 
can go to the longitudinal rest frame of the heavy system with a Lorentz boost, perform 
a transverse boost of the system such that its transverse momentum vanishes, and use 
variables that describe the kinematics of the heavy system in this frame. 



6. MiNLO merging for more complex processes 

In the present work, we have dealt with relatively simple processes, i.e. a colour-neutral 
massive particle in association with one jet. The results of ref. [14] strongly suggest that 
our procedure may be generalized to higher jet multiplicities. Since, at the level of one 
associated jet, we had to improve the MiNLO prescription with the inclusion of certain 
NNLL terms in the Sudakov form factor, it is clear that, in general, we need to find a 
similar improvement for the more general case. 

A first question that we would like to consider is whether our MiNLO procedure down- 
graded to the LO level (that is to say, to the CKKW procedure [40] as applied to the 
inclusive sample), already achieves our goal at LO accuracy, that is to say, it is such that 
by integrating the softest emission one gets a LO accurate matrix element for one less 
emission. It is easy to convince ourselves that as soon as we deal with matrix elements 
involving more than four coloured particles (including the softest emission), this is not the 
case. In fact, after the first clustering, that in our procedure simply sets the Qq scale, 
we are left with four or more coloured partons. Soft gluon resummation, in this case, 
also involve soft, non-collinear terms that arise from interference of the emission from the 
coloured external lines [41, 42, 43]. These terms are of NLL accuracy, and, according to 
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our counting, they can contribute terms of relative order y/ots. By not including them, we 
thus introduce an error of this magnitude, while LO accuracy requires the neglected terms 
to be of order a s . Notice that this problem does not manifest itself in the HJJ case, since 
precisely four coloured partons are present here (two in the initial state and two in the 
final state). In processes like ti production in association with one jet, such terms would 
have to be accounted for. On the other hand, it is possible to compute these terms using 
standard resummation techniques. It is thus conceivable that the LO MiNLO formula can 
also be improved including these interference terms. But we also stress that the CKKW 
procedure, as is, does not satisfy our requirement at leading order. 

In the present work, we have used as a clustering variable the transverse momentum 
of the boson, since it is a simple variable and the corresponding resummation formulae are 
well known. As an alternative, we could have used the hardest jet transverse momentum, 
taking B2 from ref. [44]. If this method is to be extended to processes with more than 
one radiated parton, it is clear that other clustering variables should be chosen, likely with 
good resummation properties. One should then seek either an NNLL extension of the 
CKKW procedure, or, construct a product of standard soft resummation factors, accurate 
at the NLL level, modifying the Sudakov form factor in each one to include the B2 term 
relevant to the associated configuration of clustered particles. Notice that within our 
method it was never necessary to know explicitly the hard resummation correction that is 
usually included in NLL and NNLL resummation formulae, since the NLO matrix elements 
intrinsically provide this. 

7. Conclusions 

In the present work, we have illustrated a method for constructing NLO+PS generators 
for the production of a heavy system accompanied by a radiated parton, such that, when 
integrating over the parton's phase space, one recovers the accuracy of a corresponding 
NLO+PS generator for the production of the heavy system alone. In essence, in our 
method, we start from the MiNLO prescription of ref. [14], and we look explicitly for the 
places in which it needs to be modified in order to maintain NLO accuracy on integrating 
out all radiation. We have found that the inclusion of the A2 and the B2 term in the 
resummation formula is enough to achieve this goal. Since we do not know the B2 term 
for the type of clustering that we performed in ref. [14], we have modified the prescription 
so that the transverse momentum of the boson is used instead. More specifically, the 
prescription is applicable to any clustering scheme in which the last step is performed 
using the transverse momentum of the boson as the clustering scale. 

We have not explored, at this moment, any issues related to the logarithmic accuracy 
of our approach. In other words, the logarithmic accuracy should be at the same level as 
that of the NLO+PS generators we are referring to. We postpone to future studies, a more 
accurate assessment of, and possible improvements to, the precision of the resummation. 

We have tested our method in the framework of H/W/Z production. We find that 
the method performs remarkably well. We have also found that the usual scale-variation 
method used in order to determine uncertainties may be deceiving in our case, especially for 
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processes like W/Z production, where, at the Born level, there is no renormalization scale 
dependence. We track this problem to the fact that, in the corresponding BJ calculation, 
a renormalization scale variation is already possible at the Born level. 

We point out that, using our BJ-MiNLO generators, it is actually possible to construct 
an NNLO+PS generator, simply by reweighting the transverse-momentum integral of the 
cross section to the one computed at the NNLO level. We postpone a phenomenological 
study of this method to a future publication. 

While this work was under completion, a publication has appeared that has some 
points in common with our work [11]. This work focuses on the accuracy of the matching 
conditions, and sets up a framework, using NNLL accurate resummation of soft-parton 
emission, such that one has no loss of NLO accuracy when matching. Also our approach 
is motivated by the requirement of preserving the NLO accuracy. It is more focussed, 
however, upon finding the minimal modification to the soft-parton resummation such that 
no matching is needed at all, and, in fact, we identify precisely where do we need to 
improve the resummation formula in order to achieve our goal. In this way, we find that its 
implementation in the B J case requires only a minimal modification to the MiNLO procedure, 
that, by itself, is quite simple. 
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A. The NNLL resummed differential cross section 

The NNLL differential cross section for the production of a colourless system, denoted by 

B, of virtuality Q 2 (when dealing with a single vector boson, we assume that Q is equal to 
the vector boson mass M), in the partonic scattering i + j — > B, is given by 



dg^ = a °J \- C ™ ® f«/A ( x *> j) [ C fi ® M (xb, j) x expS (c 2 Q, |) 

(A.l) 



where we have used the following definition of convolution 



[Cij ® fj/j] (x, <l) = fi/J & <?) Cij (a s (q 2 ) , |) . (A.2) 



The (To factor is such that, at leading order, 

da 



and the rapidity of the B system is given by 



I = fi/A(XA,fJ>F) fj/B(xB,fJ>F)cro, (A.3) 
/ LO 



y B = ilog^. (A.4) 
2 xb 
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The Sudakov form factor S (q 1 , q) is defined by 7 

S (q', q) = ~ [ 9 2 4 \ A («b (^)) log ^ + B («s (M 2 ))l , (A.5) 

and the Cy coefficients have the following perturbative expansion 

dj = Cij (a s , z) = 5 i3 5 (l-z) + a s c\f (z) + . . . (A.6) 

The Cij functions can thus be absorbed into a redefinition of the pdfs. We will thus assume 
now that our pdfs include this factor, and only at the end we will reinstate it. 

Equation (A.l) is easily manipulated by going to Mellin space. However, we do not 
want to lose the information on the rapidity y B of the B system. We then define the 
following Fourier-Mellin transform 

&n = f dTd^-V 2 ^-^ (A.7) 
J aq T ay B 

where 

r = xaxb , (A.8) 

so that S = Q 2 /t, where S is the hadronic center-of-mass energy. Notice that the r 
integration is performed by considering that S goes from Q 2 to infinity. We can then 
introduce a complex number N such that N = a + i (3, so that we can write 

da 



a N = C &x a &x b t Kc ^- 1 e t2lm ^y* 
Jo 

dxAdxB (xaxb) 



with 



dq%dy B 

r 1 dx A dx B N N * da 
J x A x B XA XB dq*dy B 

a J d 2 be^ rh f a/A ,N (j) fb/B,N* (j) expS (c 2 Q, |) , (A.9) 

fa/A,N (y) = J dXAX^~ 1 f a/A (xA, j) • (A.IO) 

By solving the DGLAP equations for the pdfs to the required logarithmic accuracy, and 
evolving the moments of the Cij coefficients too, we can write 

fa/A,N (y) fb/B,N* (y) = fa/A,N ( c 2<2) fb/B,N* (^Q) exp {tie \Q N (c 2 Q, y)] } 

= F N (c 2 Q) exp {Re [g N (c 2 Q, j)] } , (A.ll) 
with (see eq. (2.22) of ref. [18]) 



7 Note that this Sudakov is the square of the Sudakov form factor A 9 / 9 considered in Sec. 2 (see eq. 2.12). 
Therefore the coefficients Ai and Bi here are twice those given in eqs. (2.5) and (2.6). 
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We define 



a s = a s (clQ 2 ) , 



and, using the solution the renormalization group equation for a s , 

da s 



we write 
where 

h (y) = 
h (v) = 



-b a s - bia s - b 2 a s + ... , 
a log 

S (c 2 Q, ^)^S (a s , y) = —f (y) + h (v) + « s / 2 (y) + 

V / «q 



(A.13) 
(A.14) 

(A.15) 
(A.16) 



\ [y + log (1 - y)\ 



A 
~h 2 

M 



^log 2 (l-,) + T ^ + i^^ 
z i y 



i-y 



log(l - y) + 



y 



l - 



+ log (1 - y) , 



/ 2 (y) 



#2 y 5i6i y + log (1 - y) A 3 y 2 



bo (1 - y) b 2 



1-1/ 



(1 - J/) 1 



+ 



A 2 &! 3y 2 - 2y + (4y - 2) log (1 - y) , Ai 1 



+ 2 



&0&2 (1 - V? + b\y (1 - y)] log (1 - y) - 36 & 2 y 2 + &?y 2 + 2^ 2 y}. (A.17) 



263 



(i-y) 5 



Notice that the coefficients A 2 , B\ and B 2 do depend explicitly upon c\ and c 2 (see ref. [24]). 
The values c 2 = 1 and c\ = 2e~ lE give the coefficients that are usually reported in the 
literature. 

Similarly, we can expand Qn {c 2 Q,c±/b), using 

7aK«s) = «s7i,JV + «i72,Af , 



we can write 



Qn (c 2 Q, y ) = Qn (a s , y) = gi,N (y) + «s y 2 ,7v (y) 



with 



9i,n (y) = -^r- log (l - y) 



52,iv (y) 



&§ (i - y) 



hi,Nh log (1 - y) - (6 72,Af - &i7i,Jv) y] • 



(A.18) 
(A.19) 

(A.20) 
(A.21) 



We first perform the angular integration in eq. (A. 9), we change the integration variable 
to b = bq T and we integrate by part, to get (see ref. [18] for more details) 



daw 
M 2 



d_ ' x 



' ' (<2( ^ dg2 J ^bJi(b) expS ^c 2 Q,ciy ) exp^jy ( c 2 Q,a 



(A.22) 
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at NNLL accuracy. We then write 



y = Vo + a s b l b , 



with 



yo = -a s b log 



l b = log -J ■ 



(A.23) 
(A.24) 



Now we have to consider h as not being parametrically large, while yo is the variable that 
becomes of order 1 in the small q T limit. We thus expand S and On at NNLL order 

Sn ("S, y) = S N (a s , y ) + fo (yo) boh + y /o (yo) $1% + a s f[ (y ) b h , (A.25) 



We get 

d(TN _ 

dg2 



On (a s , y) = On ("s, yo) + «sffU (2/0) 6 Zfo • 



d f°° - 

Fat (c 2 Q) ^-exp5(a s ,y ) exp {Re [C/at (a s ,yo)] } / d6Ji(6) 
a^T jo 



(A.26) 



x exp [/ (y ) b l b + a s /( (yo) Wfc + « s Re [g' 1>N {yo)] boh + y /o (yo) ^ 6 2 J • (A.27) 
We now work out the integral 

/•OO 

/ = J dbJi(b) exp |/q (y ) Mb + a s f[ (y ) b l b + a s Re [# 1)JV (y )] hh + y /o (yo) &o^} 



1 ( £i fi ^ 

d5 Ji(5) |l + a s 6 [/( (yo) + Re y ltN (yo)] } 2^ + f / ' (y ) ^^2}^ 



1 + 2a s 6 



9 



/((yo) + Re[ 5 i,iv(yo)] ^ + 2a s / '(yo)& 2 



where 



We then obtain 



h = 2/o (yo) bo = - 



d \ f 2\ h T(l + h/2) 



] 8h 2 j \ Cl J r(l-/i/2)' 

(A.28) 



2^4i y 
&o 1 - yo ' 



(A.29) 



1 + a s fy) 



/i (yo) + Re [g' hN (yo)] 



Y>o 1 1 + J 1 + Y>o ( 1 



+ y/o' (yo) 60 



/1 



^1 1 + o + # 1 + 



+ 2^o 1 



/1 



^0 1 + 



+ 47 B Vo 1 + 



^ 1 



/1 



' r(i-|)' 



where 



^o(z) = -r(z) 



(A.30) 



(A.31) 
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Notice that the iV-dependent term, that prevents rewriting the luminosity in x space, 
depends upon N through the factor g[ N . Thus, the A^-dependent term can be factorized 
with the required accuracy as 

exp {Re [Q N (a s , y )] } exp {a s b K Re [g' ljN (y )] } 

GN{a s ,yo) a s b Kg' 1N (y ) g N *(a s ,y ) a s b K g' 1N , (y ) 
= exp exp exp exp 



ex P { \ [Gn ( a s, Vo) + a s 9i,n (vo) bo log (exp K)] 
x ex P { \ [Gn* (a s , ?/o ) + «s g[,N* (Vo) b o log (exp K)] 



exp 



\Qn\C2Q, 



K 



q T exp 



exp { -Qn* ( c 2 Q, q T exp 



K 



where 



K = Vo 1 + 77 + V'O 1 



h 



(A.32) 
(A.33) 



and we have used eqs. (A. 26) and the definition of Qn in eq. (A. 12). 

We can now perform the inverse Mellin/Fourier transform back to x space, obtaining 

d r ^ ^ , ,( K\ { „ ^ , ,( K\ 

xb, exp — J exp<S (a s , y ) J , 



d#|dy B dg^ 



[Cio ® f a / a] (xa, Qt exp y^j [C jb 



(AM) 



where Jo is equal to eq. (A. 30) with the g[ N term deleted 



Jo = 1 + a s 6o/{ (yo) 



/i 



Vo 1 + 77 + Vo 1 - 77 + 27E 



+ y/o(z/o)&o 



'■i U + ^)+V6M i 



+4 7B Vo ( 1 + -J - Vi ( 1 
+ hi-'-o ( 1 - 7j ) + ! 



2 



+ Vo 2 (l 



+ 2Vo M - 77 ) < M : 

2 



r(i + §) 



(A.35) 

r(i-|) 

At the moment we are not interested in developing this expression, that is interesting of 
its own, any further. For our purpose, we just need to expand Jo to get 



1 + 2a s 6 /( (y ) ^ + 2a s / ' (y ) b 2 



d 

1 + 2a s bof[ (yo) + 2a s /o (yo) fcjj^p 



d_ 
dh 
d 



2\ h T(l + h/2) 



a J T(i-h/2) 



12 



+ O (h 5 ) 



a s 



bof[ (yo) ( 3 h 2 - asfo (yo) b 2 ( 3 h + ... 



12 2 

Thus, the factor Jo corrects the formula in ref. [45] by terms of the following order 

3 powers of J down from leading term a s L ; 



(A.36) 



a s /i 2 « a^J 2 , 



4 powers of J down from leading term; 



2 r 4 



3 powers of J down from leading term a s L 



(A.37) 
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Furthermore, the scale choice in the pdfs, g T exp K/2, induces a change of relative order 

a s K » a s h 2 « a|L 2 . (A.38) 

If we are only interested in terms of relative order a 2 with respect to the Born contribution, 
the only term that we need to keep is the one of order a s h, i.e. the term 



-as/o (Vo) blC 3 h^ 2( 3 Ala 2 s log 



4Q 2 



This is the term kept by Ellis and Veseli in ref. [17], and the replacement 

B 2 ^B 2 + 2( 3 A 2 
is sufficient to incorporate this term. 



(A.39) 



(A.40) 



B. Scale variation in the resummed expression 



We consider the resummed cross section written in the process-independent form of eq. (13) 
of ref. [46] 

da I ,,, ,; r i / C| \ r ,.. ., , / f I \ / . a' 

b 



dq%dy B 
where 



a F j d 2 be^ rb [C ia ®f a/A ] (xa,j) [C jh ®f b/B ] ( x b,j) ex P 5(c 2 Q, 



a F = a xH, H = l + a s (Q 2 )#i + 
We further rewrite this formula as 



(B.l) 
(B.2) 



da 



dq%dy E 



a F / d 2 be^ 



C ia (p^) <8 f a / A Of) {XA 



c 



fJ-F) <8 fb/B^F) (xb) exp5 [c 2 Q, -r ) . (B.3) 



Cl 



Since eq. (B.3) is process independent, it must be possible to introduce the scale variations 
independently in the integrand and in the prefactor ap. Furthermore, in the prefactors 
involving the pdfs there is a built-in scale independence as far as the factorization scale 
variations are concerned. Thus, the Sudakov exponent must be written in a scale-invariant 
form by itself to recover the full scale dependence. We thus write 



(«*t) 



exp 



d/x 2 

2 



lcl/b 2 A* 

Performing the change of variable 



we get 



A* 



// = K n fi, 



(B.4) 



(B.5) 



exp < - 



K n c lQ 2 d/j,' 



il 



12 



A I a s 



A^ 
K 2 



log 



K 2 c 2 Q 2 



12 



+ B[a 



A^ 
K 2 



(B.6) 
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We relabel jj! — > fi, and, using 



«s {^/Kl) = a s (fi 2 ) - b log ol\ (^ 2 ) • • • , 

R, 



(B.7) 



we find 



A a s 



A ias ((j 2 ) + (A 2 + 25 ^i log K R ) a 2 (fj 2 ) 

+ [A 3 + A 1 b 2 (26i log K R + 4 log 2 K R ) + 46 A 2 logK R ] af (,u 2 ) + O (a*) , 
Bias <y ) + (B 2 + 26 Si log K R ) a 2 (/i 2 ) + O {a 3 s ) . (B.8) 



5 (^a s 
Denning 

A (q s (/i 2 ) , K R ) = A ias (fi 2 ) + {A 2 + 2b A 1 log K R ) a 2 ( / u 2 ) 

+ [A 3 + A^l (26! log K R + 4 log 2 if R ) + Ab A 2 log tf R ] a| (/x 2 ) + O (af ) , 
5 (a s <y) , if R ) = Bia s (/i 2 ) + (B 2 + 26 £i logtf R ) a 2 ( M 2 ) 

+ [Aia s (fi 2 ) + (A 2 + 26 Ai log K R ) a 2 s (/x 2 )] 21ogK R 
= [Bi + 2A 1 logif R ]a s (fi 2 ) 

+ [B 2 + 2 (A 2 + b B 1 ) logK R + 46 v4i log 2 K R ] a 2 (^ 2 ) + O (af) , (B.9) 



we get 



exp 



A (a s ((m 2 ) , K R )) log + B (a s (^ 2 ) , K R ) 



We can now break the integral in the exponent into two integrals 

^ 2 c 2 Q 2 ,c 2 Q 2 /-^Q 2 

/*T R c 2 /& 2 "' Jk14/^'" J<$Q> 
and we call the last integral / 



/' 

JK 



I 



d^ 2 



C 2 Q 2 ^ 



A (a s ((/i 2 ) , K R )) log + B (a s (// 2 ) , K R ) 

A 1 



(B.10) 
(B.ll) 

(B.12) 



Since this integral does not have large logarithms, we can evaluate it at order a s (c^Q 2 ), 
ignoring the running of a s , 



= Aia s (c 2 Q 2 ) ^ log ^ + (Bi + 2Ai log K R ) a s (c 2 Q 2 ) / 



clQ 2 d/U 2 



/^c|Q 2 A* 

= a s (clQ 2 ) [2Ai log 2 K R - 2 (Bi + 2Ai log K R ) log K R ] 
= -2 a s (c^Q 2 ) [Ai log 2 K R + Bi log K R ] . 

Using the renormalization group equation for a s 

dfi 2 _ 1 



^ rC |Q 2 A^ 



2 dog 



a 



S 



da s , 



(B.13) 



(B.14) 
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we can write the identity 



^da s (/i 2 ) = — a s (clQ 2 ) - a s 



b 2 



. (B.15) 



Solving this equation for a s (c^Q 2 ) and inserting it into the expression of / of eq. (B.13) 
we get 



L 



- as 



Kef 

b 2 



a 2 s (fi 2 ) [-2b A 1 log 2 K n - 2b B 1 log K B 
[2A 1 log 2 K R + 2B 1 logK R ] . 



(B.16) 



Using eqs. (B.ll), (B.12) and (B.16), we can write eq. (B.10) as 



S 



(<*• 7) 



cxp ■ 



/; 



V L t_ u.2\ T, u, e\Q 2 



x <! 1 - a s 



A (a s ((fi 2 ) , K R )) log - 2 -| — I- B' {a s (fi 2 ) , K R ) 



K 2 r 2 

R 1 

b 2 



[2Ai log 2 K R + 2Bi logK R ]|, 



(B.17) 



where 



5' (a s (^ 2 ) , K R ) = [Bi + 2A 1 log K R ] a s (/i 2 ) + [B 2 + 2.4 2 log K R + 2 fe A 1 log 2 if R ] a 2 (^ 2 ) . 

(B.18) 

The last term in curly braces of eq. (B.17) affects the Cij terms by a contact term. In our 
case it is irrelevant, since we always deal with the O (a 2 ) expanded cross section multiplied 
by the Sudakov exponential. Notice also that the argument of a s in the contact term can 
be changed by a factor of order 1 at the required accuracy. Thus, the K R dependence drops 
from there, and a s can be evaluated at a scale of the order of the low scale, i.e. c\/b 2 . 

B.l Process dependent form 

We can also translate eq. (B.18) to the case of a process-dependent form, which is similar 
to formula (B.3) except that the prefactor <rp is replaced by do. In order to perform this 
translation, the H factor must be somehow incorporated in the Cij coefficients and in the 
Sudakov exponent. In order to do this, we first notice that if the Born term is of order a™ 
in the strong coupling constant, the scale dependence of H is derived using the identity 

off (i^ 2 Q 2 ) H F C (a s (K 2 R Q 2 ) , K R ) = < (Q 2 ) H F C (a s (Q 2 ) , l) , (B.19) 
which, using eq. (B.2), leads to 

H F (a s (KlQ 2 ) ,K R )=l + a s (kIQ 2 ) [Hi + nb log (B.20) 
This can be replaced by 



H F (a s {KlQ 2 ),K B ) 



1 + a s 



Kle 2 
b 2 



{Hr+nb^ogKl) 



x exp <^ (Hi + nbo log K R ) 



a s {KlQ 2 ) - a s 



K 2 r 2 

A R C l 

b 2 



-,(B.21) 
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and using eq. (B.15) 



I ; r, s ( K ffi ) (H 1 + 2>nb Q ]ogK VL) 



6 2 



x exp 



{ - / -V«i (/u 2 ) 6 [ffi + 2n6 log K R ] } . (B.22) 



l/b 2 V 



On the other hand, the scale dependence in the resummation factor must be equal to what 
we have found in eq. (B.17). Inserting eq. (B.22) in the process independent formula, we 
should get the process dependent result. This induces the modification 8 



B% = B 2 + 6 #i + 2nb 2 log K R , 



and the modification of the Cy coefficients by a contact term 

dj (z) -> dj (z) + S(l-z)^(H 1 + 2nb log K R ) a s ( 



KM 
b 2 



(B.23) 



(B.24) 



where the 1/2 factor comes from the fact that, in the resummed expressions, there is the 
product of two Cij terms. 



C. A mathematical complement 

In this section, we explicitly estimate the size of the following integral, that we use through- 
out to estimate the contributions to the inclusive cross section 



with 



«s(m 2 



bo log 



where A is the usual A QCD . We first evaluate the argument of the exponent 



j _/ 



v a 2 ; / ;n i () g" 

A <- L 
bo 



2 

A2 



where we have defined 



I = lOB 



A 2 



L = log 



91 

A 2 ' 



x = log 



A 2 ' 



(C.2) 



= U l L < x ^ L - x)= £( L]o *i- L+i )' (c - 3) 



(C.4) 



8 Note the difference of a factor 2 between the last term in eq. (B.23) and the corresponding factor in 
eq. (3.8), due to the square factor in the definition of the Sudakov form factor we are using in eq. (2.12). 
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Equation (C.l) then becomes 

-L 



I(m, n) 



f dl(L- l) m exp {-a [L (log L - log I) - L + /]} 
Jo °o 1 

-I' 
-[' 



dl exp{mlog(L — I) — nlog/ — a [L (logL — log/) — L + I]} 



where 



dl exp [/(/)] , 



/(/) =mlog(L-/)-nlog/-a[L(logL-logZ)-L + Z] , 



(C.5) 



(C.6) 



and we have defined a = A/ bo. This integral has to be computed for large L. We look 
for an approximation of the integral using the saddle-point technique, i.e. we expand the 
argument of the exponent around its maximum 



f(l) = f(ht) + \ f"(l M ) (I - Im? + 0({l- l M f) 



with 



Im 



2a 



m — n 



+ 2aL — \fjrn — n) 2 + 4 am L 



= L- 



mL 



+ 0(1) 



For large L we have 



f(l M ) = (y-n) logL + 0(1) , 
f"{l M ) = ~ + 0(L-^) , 



(C.7) 

(C.8) 

(C.9) 
(CIO) 



and we get 



I(m,n)^ -Lf- n d/exp --{l-l M ) 2 = L? ~ n / d/'exp 
°o Jo 1 L J °o J-i M 



L 



L (L-l M ) y \ m+1 r\f™ 

dxexp(— x 2 ) m——=L^ n / dxexp(— x 2 ) 

J —\faL 



b%y/E 



< 



m+1 

2 . (C.ll) 



In other words, for each power of the logarithm in eq. (C.l), we lose half power of a s . 
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